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SUMMARY 


A study is made for the development and computations of the separating 
and reattacing shear flows. The highlight of the study is an attempt to 
predict the third-moments of turbulent velocity which is responsible for the 
diffusion transport of the Reynolds stresses. The present computations show 
that the third-moments obtained by employing the low-Reynolds number model of 
transport equations improve the prediction of the third-moments. 

The modeling for scalar variables are also performed for the heat transfer 
computations. Since the transport equations for <u.e> and <u.u.e> have been 

» J 

given in the last NASA CR, the study has been extended further to the modeling 

2 2 2 
of <e >, <u.e > and c Q (dissipation rate for <e >). The formulations are 
1 0 

shown in this report. 
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INTRODUCTION 

Separation and recirculation of flow is encountered in a vast array of 
engineering applications. This phenomenon of flow separation and 
recirculation associates itself with higher turbulence levels which not only 
render the flow greater analytical complexity but also would result in highly 
augmenting its heat transfer and momentum aspects. Thus, demand for 
mathematical models that can predict any complex turbulent flows is increasing 
as well as experimental measurements for better understanding of such complex 
flows. It has been observed by many researchers [1,2] that non-isotropic 
turbulence effects are quite significant in separating shear flows and 
consequently occurring reattaching and recirculating flows beyond a step. 
Moreover, turbulent third-moments and corresponding triple scalar-velocity 
products are equally predominant in the separated shear flow field interacting 
with turbulent boundary layers. 

The most noticeable features of the experimental observations reported by 
Chandrsuda and Bradshaw [1] and later by Driver and Seegmiller [3] are that 
the shapes of the third-order moments of turbulence fluctuating velocity 
change rapidly along the separated shear layer due to distortion of the large 
eddies by effectively irrotational mechanisms. The implication of these 
observations in separating and reattaching shear layer is that algebraic 
relations between the third-order moments and the Reynolds stresses are likely 
to be inadequate for strongly perturbed flows, making it necessary to use 
transport equations for the third-moments. 

This study is conducted to develop new models of turbulence with the 
third-order closure in order to predict the turbulent flow through a duct with 
a relatively complex geometry. As a basic study of modeling, a backward-facing 
step is chosen to carry out computations. The advantage of having such a flow 
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geometry is that a number of experimental data are available in the literature 
for comparison with present results. Computations of the momentum and 
temperature fields in the flow domain being considered entail the solution of 
the time-averaged transport equations containing the second-order turbulent 
fluctuating products. The third-order products, which are responsible for the 
diffusive transport of the second-order products, attain greater significance 
in separating and reattaching flows. The following chapter describes the 
modeling of these transport equations. 


MATHEMATICAL MODELS FOR THIRD-MOMENTS 

The transport equations for the Reynolds stresses can be given as 
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where P. . and <|>.., respectively, represent the production and the 

J J 

pressure-strain rates. In the parenthesis of the above equations, the upper 
case letter U and the lower case letter u represent the Reynolds averaged and 
the instantaneous values of the velocity, respectively. The symbol < > 
denotes that the argument inside this bracket is ensemble-averaged. The 
pressure-strain correlation <t>. . is evaluated by using the model of Launder 
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et al. [4]. 

The diffusion rate, D.., contains the terms given as follows: 
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where p, p, and v represent, respectively, the turbulent fluctuating 
pressure, the fluid density and the fluid kinematic viscosity. The symbol 6 
denotes the Kronecker delta. Since the diffusion rate of the Reynolds 
stresses is governed by the gradient of the third-moments <u ^ u j u | ( > * it 
becomes necessary to evaluate the third-moments exactly. Thus, the 
transport equations for such third-moments are developed and shown as follows. 

The third moments, <u i u j u | < > * can be evaluated by formulating 
their transport equations as 
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In the above equations, terms I and II represent, respectively, the 
generation rates due to mean strain and the Reynolds stresses. Term III is 
responsible for the generation rate if the Gaussian approximation is applied; 
this can be merged into term II. Term IV, which is the pressure-stress 
correlation, plays a significant role in redistributing the tensor and in 
diffusing due to fluctuating pressure. Term V can be divided into the 
molecular diffusion rate and the dissipation rate. 

In order to account for the viscous effects that are predominant near the 
wall, a low-Reynolds number modification needs to be incorporated in the 
transport equations of the third-moments. 

The theory behind this low-Reynolds number model is based upon the fact 
that the scale of the eddies created in the shear layer changes rapidly as the 
flow approaches the solid wall by enhancing the turbulence energy dissipation 
rate. Thus, the model should be devised to control the size of the eddy. 

This was done by evaluating the correct expression for the dissipation rate, 
c, by formulating it in each layer near the wall such as the viscous 
sublayer, the buffer layer, and the fully turbulent core flow. This 
controller is incorporated in the pressure-stress term (Term IV in Eq. (3)) in 
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order to adjust the time of rate change of the third-moments of turbulence 
velocity as follows: 


6 ijk ~ Pressure-stress term 


- C T k <U l“jV 


and 
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e... .. = near-wall correction for pressure-stress term 
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where k and y represent the turbulence energy and the distance from the wall, 
respectively. 

The dissipation rate is given as: 


c ijk s Dissipation due to viscous effect 
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= Diffusion due to viscous effect 
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The values for the constants are listed in Table 1. 
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TABLE 1. Recommended values for the constants 


C 9 

C 4 


c 

Y 

C Y w 

C 

cy 

1.0 

2.55 

3.0 

5.8 

( low-Re no. model) 
(high-Re no. model) 

8.0 

0.2 

The models 

considered 

in this 

study are defined as 

follows: 



High-Revnolds number model 
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P ijk 2 = Pr °duction due to Reynolds stresses interacting with their 
gradients 
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Low-Revnolds number model 
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SOLUTION PROCEDURE 


Figure 1 shows the solution domain of the flow field. The inlet section 
is located approximately two step heights upstream from the step and the 
outlet section about sixty step heights downstream therefrom. After several 
grid tests it was discerned that the system with 62 x 62 node points gives 
results independent of the grid. More details are given in the following 
chapter. The computations were performed using the 62 x 62 variable grid with 
the grid expanding linearly at the rate of 2% in x-direction and at 3% in 
y-direction. Figure 2 shows the numerical grid used in computations (5:1 
expansion in y-direction). 

Computations of the transport equations described in the preceding 
section are achieved by using the finite volume method [5]. The iterative 
procedure is terminated when the maximum value of the relative residual 
sources of U, V, and mass balance falls below 1%. However, the computations 
of the triple products are continued until the relative residual sources fall 
below 3 x 10 -8 . 

RESULTS AND DISCUSSION 

Figures 3 and 4 show the computed skin friction coefficient, C^., and 
pressure coefficient, C p , respectively, along the bottom wall. The results 
are also compared with the experimental data of Driver and Seegmiller [3]. 
These coefficients are defined as follows: 
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where t w and P represent the wall shear stress and the Reynolds averaged 
pressure, respectively. The subscript IN stands for inlet and ref indicates 
the reference location which is two step heights upstream of the step. 

Computations have been performed with several different grid systems 
ranging from 32 x 32 to 62 x 62. It is shown in these figures that fair grid 
independent status is attained with the grid system finer than 52 x 52. For 
this reason the grid with 62 x 62 has been used in the rest of the 
computations. The agreement between the computation and the measured data is 
reasonably well except the discrepancy for the pressure coefficient beyond the 
reattachment point. 

Figure 5 shows the mean velocity profiles at two downstream locations. 

The computed profiles are compared with the experimental data taken by 
Chandrsuda and Bradshaw [1]. The measured data were obtained by using two 
different methods: one by pressure-probe and the other by hot-wire. 

Agreement between the computation and the two sets of the experiment is quite 
reasonable. It is also noticed that the trend of the computed results accords 
more with the hot-wire measurements than with the pressure-probe data. 

Figure 6 shows the profiles of the Reynolds stresses at several locations 
downstream of the step. As seen in this figure agreement of the normal 
stresses (<uu> and <vv>) and the shear stress (<uv>) between computed results 
and the measurement is fairly well within the error of 10% beyond the 
reattachment position; however, the computed <uu> and <uv> profiles within the 
recirculating region do not agree with the measured levels at the peak 
positions. This disagreement seems to be larger near the step, but both 
results coincide toward far downstream. 

Probably, the generation rate and the pressure-strain correlation 
associated with the secondary mean strain rates (av/ax) for the transport 
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equation of <uv> is not as high as it should be in this near step region, 
whereas the actual shear stress level might be enhanced due to the cornering 
flow at the bottom of the step. 

Figure 7 represents the comparison of the third-moments of turbulence 
fluctuating velocity between the high-Reynolds number and low-Reynolds number 
models. In this figure large changes are observed between these two models. 
The high-Reynolds number model gives higher levels in the near-wall region 
than the low-Reynolds number model. This fact indicates that the terms 
represented by Eqs. (5) and (6) play a significant role in the determination 
of the third-moments. In the near-wall region, the eddies created through the 
separated shear layer tend to become smaller, and as a result, the rate of 
turbulence energy dissipation increases rapidly. This is why the effect of 
c added in the model promotes the dissipation action of <u.jUjU k > in 
the near-wall region giving lower levels of the third-moments toward the 
solid wall. Agreement of the computations by the low-Reynolds number model, 
thus, is considerably improved. 

Figure 8 represents the comparisons of the low-Reynolds number model with 
other algebraic models developed by several researchers: Hanjalic and 

Launder [6], Cormack et al. [7], Daly and Harlow [8], and Shir [9]. As shown 
in this figure, the results obtained by the transport equations for the 
third-moments give generally better agreement with the experimental data than 
the results obtained by using the corresponding algebraic models. However, it 
is also noticed that the levels of the third-moments predicted by using the 
transport equations are very high across the separated shear layer above the 
recirculating flow, although they agree reasonably well with the experimental 
data at the location near the reattachment and the downstream region. 
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Another advantage of using the transport equations is that a low- 
Reynolds number correlation can be incorporated in order to account for the 
viscous effect near the wall as well as the convection effect. Thus, the 
prediction by using the transport equations is closer to the measured data 
near the wall throughout the flow region. 

The above discussion can be summarized as follows: 

1. The Reynolds-stress model developed in this study is shown to predict the 
separating and reattaching shear flows properly. 

2 . The third-order turbulence velocity fluctuating tensor denotes rapid 
changes in the reattaching and recirculating flow regions. This behavior 
can be predicted by developing the transport equations for the 
third-moments of turbulence velocity. 

3. It was shown that by promoting the dissipation effect of the third- 
moments in the near-wall region the prediction of the third-moments were 
considerably improved. This low-Reynolds number model for the 
third-moments provides more universal results than do the algebraic 
models. 

Figures 9-11 show variations of the computed Reynolds stresses. Figures 
12-15 show variations of the triple products of the velocity fluctuations. It 
is clear that the sharp changes in the levels of the triple products are 
produced in the separating shear flow region which diminish toward downstream. 

DEVELOPMENT OF TEMPERATURE FLUCTUATING EQUATIONS 

A new model has been under developing for the computation of heat 
transfer rates in a separating and recirculating flow. Since a transport 
equation approach is required for better prediction of such complex flows, 
temperature fluctuating variables with interaction with velocity fluctuating 
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transport are correlated with defined hydrodynamic and Reynolds averaged 
velocity fluctuations. 

In the last report [10] the formulations and the preliminary computed 
results were shown for the variables of <u.e> and <u^u^e>. In this 
report formulations of other variables such as 

R ee = <e 2 >/2 (14) 


R iee = <u i e > 


(15) 


and 


,36 36 . 

c = (X< ? 

6 3X.3X. 

vl J 


( 16 ) 


are presented where a and 6 denote, respectively, the thermal diffusivity 
and fluctuating component of temperature. c_, which represents the 
dissipation rate of turbulence temperature energy level, corresponds to e. 

The transport equation for R.. was developed because it appears in 

UU 

the buoyant source term of the vertical heat flux equation of stratified 
flows. For this reason the transport equations for c a and R. were 
developed in order to close the equation. The buoyant term is 
important because a fluid element receives an upthrust (buoyant effect) 
proportional to a6g if its instantaneous temperature exceeds the average 
temperature at the horizontal plane by e. This will produce an extra 
enthalpy flux proportional to a6g(T + 6) which upon time averaging 
reduces to agR_ and increases the overall levels of <ve> [11]. 
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Scalar Double Product R 
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Transport equation for R is given as: 
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Equation (17) is explained as follows: 

(i) Rate of change of turbulence temperature energy level 
balances with 

(ii) Production of turbulence temperatures energy by mean temperature 
gradient 
minus 

(iii) Energy dissipation rate of turbulence temperature by molecular action 
plus 

(iv) Diffusion transport due to turbulence velocity fluctuation 
plus 

(v) Diffusion transport due to molecular action 

The dissipation rate for the double product of the temperature 
fluctuations, c 0 , is obtained in a similar manner as is done for c. 
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where 


i 



a 


ae 



ae 



Unaveraged quantity 


Equation (18) is interpreted as follows: 

(i) Rate of change of c a 
balances with 

(ii) Generation due to mean temperature gradient 
plus 

(iii) Secondary generation due to mean temperature gradient 
plus 

(iv) Generation due to mean vortex stretching action 
plus 

(v) Generation due to turbulence fluctuations 
minus 

(vi) Destruction of c Q due to fine scale turbulence interactions 

y 

plus 

(vii) Diffusion rate due to turbulence velocity 
plus 

(viii) Diffusion rate due to molecular action. 

It is generally recognized that the contribution of terms (ii), (iii) 
and (iv) is negligibly small in high turbulent flows. Thus these terms can be 
dropped. 

Terms (v) and (vi) are treated in the same manner as c: 
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where 


P = - <u. e> — 
e k ax. 


( 20 ) 


Term (vii) may be given as 
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Equation (23) may be explained as 

2 

(i) Rate of change of <u^e > 
is balanced with 

(ii) Generation due to mean temperature gradient 
plus 

(iii) Generation due to mean vortex sketching 
plus 

(iv) Generation due to turbulence 
plus 

(v) Diffusion due to turbulence 
plus 

(vi) Diffusion due to thermal molecular action 
plus 

(vii) Diffusion due to momentum molecular action 
plus 

(viii) Pressure-scalar correlation 
minus 

(ix) Buoyancy effect. 

The pressure-scalar correlation is split into two terms as 
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(24) 


The buoyancy term is given as 

(ix) = - <e 2 > g s. 2 (25) 

o 

Thus, by neglecting molecular diffusion terms we can solve Eq. (23) after 
determining coefficients appearing in diffusion (term (v)) and pressure-scalar 
correlation (Eq. (24)). 
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FIGURE CAPTIONS 


FIGURE 

1 . 

Solution domain — Backward-facing step geometry with the 
separating, recirculating and redeveloping regions. 

FIGURE 

2. 

Numerical grid used in computations. 

FIGURE 

3. 

Skin friction coefficient along the bottom wall downstream 
of the step. 

FIGURE 

4. 

Pressure coefficient along the bottom wall. 

FIGURE 

5. 

U-velocity profiles downstream of the step. 

FIGURE 

6. 

Reynolds-stress profiles downstream of the step. 

FIGURE 

7. 

Third-moment profiles downstream of the step: Comparison 

of the low and high-Reynolds number models. 

FIGURE 

8. 

Third-moment profiles downstream of the step: Comparison 

of the low-Reynolds number model with algebraic models. 

FIGURE 

9. 

Variation of <u^> 

FIGURE 

10. 

Variation of <v?> 

FIGURE 

1 1 . 

Variation of <uv> 

FIGURE 

12. 

Variation of <u^> 

FIGURE 

13. 

Variation of <u^v> 

FIGURE 

14. 

Variation of <uv 2 > 

FIGURE 

15. 

Variation of <v^> 
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FIGURE 1. Solution domain— Backward-facing step geometry with the 
separating, recirculating and redeveloping regions. 
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FIGURE 4. Pressure coefficient along the bottom wall. 
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FIGURE 8 
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FIGURE 9. Variation of <u 2 > 
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FIGURE 10. Variation of <v 2 > 
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FIGURE 11. Variation of <uv> 
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FIGURE 12. Variation of <u 3 > 
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FIGURE 13. Variation of <u 2 v> 
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FIGURE 14. 


Variation of <uv^> 





